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Abstract 

This paper revisits Noether's theorem on the constants of motion for 
Lagrangian mechanical systems in the ODE case, with new results on 
both the theoretical and the applied side. We make full use of invariance 
up to a divergence, or, as we call it here, Bessel-Hagen (BH) invariance. 
By recognizing that the Bessel-Hagen (BH) function need not be a total 
time derivative, we can deduce nonstandard constants of motion, involving 
for example the Hamiltonian action. We prove that we can always triv- 
ialize either the time change or the BH-function, so that, in particular, 
BH-invariance turns out not to be more general than Noether's original 
invariance. We also propose a version of time change that simplifies some 
key formulas. New applications include dissipative systems, homogeneous 
potentials and new superintegrable systems. Among the more standard 
examples, we give a simplified derivation of Laplace-Runge-Lenz vector 
for Kepler's problem. 

Keywords: Noether, invariance, constants of motion, Bessel-Hagen. 
AMS subject classification: 34C14. 

1 Introduction 

Emmy Noether's theorem on conserved quantities in mechanical systems is 
widely celebrated as a fundamental tool in modern theoretical physics. Noe- 
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ther's original 1918 work [13] was motivated by general relativity, and to this 
day most applications of the theorem are in field theory. Mathematically ori- 
ented treaties that prove Noether's theorem, such as the ones by Olver [T3] and 
Giaquinta and Hildebrandt [5] , are in a partial differential equation setting and 
use advanced tools, such as generalized vector fields. Seminal papers in the field 
are the ones by Trautman [T7] and Rund [15 . For historical accounts of how 
Noether's theorem has been used, the reader can refer for example to Sarlet 
and Cantrijn's review [TB], and the books by Giaquinta and Hildebrandt [5], 
Olver [T3], Neuenschwander [T3] and Kosmann-Schwarzbach [5]. 

It is lamentable that when we come to discrete classical mechanics, most 
textbooks give short shrift to Noether's theorem. For example, Arnold [I] only 
proves an easier special case, that is not powerful enough to yield the conser- 
vation of energy for autonomous systems. There have been some authors that 
tried to restore Noether's theorem to its rightful power in classical mechanics, 
as for example Levy-Leblond jTU] and Desloge and Karch [7] . The present work 
goes further along that way, with a thorough rethinking of the theory, and a sys- 
tematization and expansion of the applications, within the ODE setting. Our 
independent variable will always be time t 6 tt, our state space will be the 
Euclidean space R™, and we will not use any geometry of manifolds. 

To start with, we have isolated in Lemma[T]a step that is in common in all the 
various formulations of Noether's theorem: it is a basic formula about how the 
Hamiltonian action changes when the first variation is allowed to have variable 
cndpoints. The simplest version of Noether's theorem (Subsection 12 .2[) will be 
an immediate corollary of that, assuming that the action is (infinitesimally) 
invariant under space change. 

Invariance "up to a divergence" was Bessel-Hagen's 1921 generalization [2] 
of the strict invariance of the original paper by Noether. Since the "divergence" 
only makes sense in a PDE setting, we are going to call this concept "BH- 
invariance" , or "invariance up to a BH-function" . For its ease of use and full 
power, BH-invariance under space change is our favourite setting for Noether's 
theorem, of which we will give three variants in Section [3J The version in 
Theorem [5] is the most natural to justify and to prove. However, we came to 
discover that it was more general than we had initially expected: it allows for 
a constant of motion that is not a point function of time, position and velocity, 
since it involves an integral. Such constant of motion looks like this: 



and does not require special assumptions on the variation q e . The formal state- 
ment is in Theorem [3] We found two neat applications of this theorem: some 
dissipative system f Subsection 13.21 and Section [6]) and the conservative systems 
with homogeneous potential (Section [7]) exhibit a constant of motion that in- 
volves energy and Hamiltonian action. These results seem to be new. 

Theorem2]is our third version of Noether's theorem for BH-invariance under 
space change. It imposes restrictions on the BH-function (the "total derivative 



t^dgL(t,q(t),q(t)) -d e q e (t) 
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L(s,q e {s),q £ (s)) _ ds , (1) 
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condition" ) and on the space change, which ensure that the constant of motion 
will be an old familiar function of t, q, q. We will follow the unifying scheme of 
Theorem 2] in most of the later examples. 

In Section Q] we show how time change can be added to Noether's theorem 
toolbox, first with Emmy Noether's original assumption of invariance with space 
and time change (Theorem O and then with BH-invariance under space and 
time change (Theorem [6]) . On the face of it, Bessel-Hagen's formulation would 
appear the most general. Boyer [5] in 1967 noticed that we can always "trivialize 
time change" with a simple transformation, and obtain the same constants of 
motion using space change and a modified BH-function. Olver Q31 Exercise 5.33] 
remarks without details that for each conservation law of the Euler-Lagrange 
equation there is a corresponding strict (i.e., with time change and null BH- 
function) variational symmetry which gives rise to it via Noether's theorem. 
Our contribution here (Theorem [7]) is an explicit formula that "trivializes the 
BH-function" if we are given space change, time change and BH-function, and 
get the results via a modified time change and null BH-function. In other words, 
Emmy Noether's original formulation is not less general than Bessel-Hagen's 
version because of a direct, constructive transformation, and not just through 
the non-elementary converse of Noether's theorem. In each single example it will 
be a matter of elegance which one to use: either time change or BH-function. In 
the system of Section [TUJ using time change and BH-function together appears 
to be best. In Section[7]we discover an example of BH-invariance using the total 
derivative condition, but then the alternative time change will turn out to have 
a nicer formula and be superior in retrospect. 

To be precise, the time change that we introduce in Subsection 14. 1[ and use 
throughout the examples, is not the one that is customary in the literature. 
In Subsection 14.31 we show that our new approach is equivalent to the tradi- 
tional one, but it simplifies some key statements, particularly Theorem [7] We 
have decided that this was reason enough to standardize our paper to the new 
definition. 

In the second half of the paper we will review several examples of applications 
of Noether's theorem to classical mechanics systems. In Section[S]we do not add 
anything new to the classical discussion of momentum and angular momentum, 
which need neither time change nor BH-function. With conservation of energy 
we highlight how the time change and the BH-function approach relate to each 
other, and also how Theorem [2] allows for a BH-function that is not a point 
function of t, q, q. 

In Section[5]we expand on the already mentioned dissipative systems. When 
the potential is homogeneous of degree 2 the constant of motion {l} becomes 
an ordinary first integral, providing a new approach to a well-known result by 
Logan [TTj . 

In Scction[7]we show that for Lagrangians of the form L = ^m\\q\\ 2 — U(q) 
with homogeneous potential (U(Xq) = X q U(q)) there is a constant of motion 
that involves energy and Hamiltonian action. In the special case k = —2 this 
constant of motion reduces to a point-function of t, q, q (formula (f54|) ). which 
has illuminating consequences for the dynamics (formula (|56[) in particular). As 
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far as we are aware, these results were well-known to Celestial Mechanics only 
in the central potential case U(q) = — fc/||g|| 2 , which is further distinguished by 
its trajectories being Cotes' spirals. 

In Section [5] we deduce the Laplace-Runge-Lenz vector constant for Kepler's 
problem. Our space variation q £ is simpler than usual in the literature, thanks 
to the use of a delay term q(t + e) in the space change. We also take the 
opportunity to show how our formulation of Noether's theorem can compute 
functions that are constant along single motions, and not necessarily along all 
solutions to Lagrange equations. 

In Section Owe use the technique of Theorem 2] to extend the class of super- 
integrable system that was found with a different method in a recent paper . 
Finally, in Scction fTOl we study the system of a particle in a plane- wavelike field, 
that we have taken from a paper [4 by Bobillo-Ares, and which is notable be- 
cause the version of Noether's theorem that can be most naturally applied to it 
is the one with time change and BH-function combined. 

Notations 

Throughout this work, R" will be the usual Euclidean space, and we will use 
the notation x ■ y for the inner product and ||x|| for the norm. The symbol q 
will denote either a vector in R n or an R"-valued function of une variable, and 
the context should clarify which is the case. Given two points x = (xi, x 2 ), y = 
(2/i)2/2) in the plane R 2 we set det(x, y) :— xiy 2 — x 2 yi- For partial derivatives 
of expressions with respect to a variable x we use the fractional notation d/dx, 
whilst for the partial derivatives of a named function F we write more simply 
d x F(x,y) or d 2 F(x,y). The gradient of a smooth scalar function q i-» f(q) of 
n variables will be denoted by d q f(q), and it will be treated as a vector in R™. 
To clarify our usage without a formal definition, here is Taylor's formula for a 
scalar function of (p, q) eM" xR": 



f(p + h,q + k) = f(p, q) + d p f(p, q)-h + d q f(p, q) ■ k + o(\\p\\ + \\q\\) . 

2 Invariance under space change only 



The variational approach to classical mechanics starts with a smooth Lagrangian 
function Lit, q, q) defined in an open subset of R x R™ x R™ with values in R, and 
trajectories q(t) such that L(t, q(t), q(t)) is defined, so that we can introduce the 
Hamiltonian action functional as 



We then posit that the mechanical motions will be fixed-endpoint-stationary for 
the action functional, in the sense of the Calculus of Variations. As well known, 
the stationarity condition is equivalent to the Lagrange, or Euler-Lagrange, 
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differential equations: 

J f dgL(t, q(t),q(t)) - d q L(t, q(t),q(t)) = Vt . (3) 

Noether's contribution was to show that mechanical conservation laws can be 
deduced from some kind of invariance property of the action functional. Let us 
review these invariances, starting from the simplest case. 



2.1 Space change and Hamiltonian action 

Definition 1. Given a trajectory q(t), we will call space change a smooth one- 
parameter family of trajectories (e,t) H> q £ (t), with values in K™, defined for 
the same times t as the motion q(t), and for e close to 0, such that q £ (t) reduces 
to q(t) when e = 0. 

This notion of space change is general enough to encompass what we call 
the time-shift family q £ (t) := q(t + e) as a special case. We still call it "space" 
change because the Lagrangian L(t, q £ , q £ ) depends on e only through the space- 
velocity slots. In Section [4] we will let the time slot depend on e too, by means 
of what we will call "time change" . 

Consider the action functional along the family q e as a function of e: 

e i ^ A a<b (q £ ) . (4) 

One crucial and easy fact is that the derivative of the action with respect to e 
at e = has an integral-free formula if we assume that q(t) is a solution to the 
Lagrange equations: 

Lemma 1 (Integral- free formula for the derivative of the action). If q(t) is a 
solution to the Lagrange equations then 

dA a , b (q £ 



r=0 



d 4 L(t,q(t),q(t)) -d £ q £ (t) 



£= 0,t=b 



-d 4 L(t,q(t),q(t))-d e q e (t) 



(5) 



Proof. Using the Lagrange equation 
d_ 



^-L{t,q e {t)AS)) = (6) 



= [d q L{t,q(t),q(t))) ■d s q e (t)\ e=0 +a i L(t,q(t),m) ' ( 9 ^eW| £=0 ) = 
= (j t d 4 L(t,q(t),q(t))) ■ d £ q £ (t)\ e=0 + 

+ d 4 L(t,q(t),q(t)) ■ J t (d £ q £ (t)\ 5=0 ) = 

= j t (d q L(t,q(t),q(t))-d e q e (t)\ e=0 ). 

Integrating with respect to t and taking the derivative g~ out of the integral 
sign we obtain equation §5§ . □ 
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If q e (a) and q e (b) do not depend on e, formula (J5J) is zero. This fact should 
not surprise, since q(t) by definition makes the action stationary with respect 
to variations that keep the endpoints fixed. Lemma Q] simply generalizes to the 
case when the endpoints of the variation are not fixed. 

2.2 Simplest version of Noether's theorem 

Classical mechanics textbooks that have chosen to dwell on Noether's theorem 
as little as possible usually present the following simple version, which only uses 
invariance under space change (Arnold's book [I], for example). 

Definition 2. We will say that there is infinitesimal invariance under the space 
change q e if 

=0 Va,6, (7) 

oe e=0 

or, equivalently, that the Lagrangian is invariant: 

|-L(t,<fe(*),fc(t)) =0 Vt. (8) 
Oe e=0 

Theorem 1 (Noether's theorem with space change only). Suppose that the 
Hamiltonian action A is infinitesimally invariant under the space change q e , 
and that q(t) = qo(t) is a solution to Lagrange equation. Then the function 

t^dgL(t,q(t),q(t)) ■d e q e (t)\ e=0 - (9) 

is constant in t. 
Proof. Let us set 

F(t) =dqL(t,q(t),q(t)) -deqe(t)\ e=0 . (10) 

From Lemma Q] we have that 



dA a ,b(qe) 



de 



e=0 



F(b) - F(a) 



Combining this with assumption (O, we obtain that F(b) — F(a) = for all a, b, 
which means that F is constant. □ 

There are situations when not only the function e ^ A a ^(qe) has zero deriva- 
tive at e = 0, but it is actually constant, or, equivalently, L{t,q e (t),q e {t)) does 
not depend on e. These will be called finite invariances under space change. 
The standard examples are the Lagrangians that are invariant under either 
translations or rotations (Examples [5] and [3] in Section [SJ, whereby any smooth 
trajectory q(t) can be embedded in a translated or rotated family q £ (t) with 
the property of finite invariance. When q(t) also solves the Lagrange equations, 
the momentum or angular momentum will be conserved respectively. Although 
very simple and neat, this concept of finite invariance does not seem powerful 
enough to cover conservation of energy, for instance. 
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3 BH-invariance under space change 



Definition 3. We will call BH-function a smooth real function G(e, t), defined 
for the same times t as the motion q(t), and for e close to 0. The trivial BH- 
function will simply be the constant 0. 

We can modify the function (U) by introducing a BH-function G(e,t), and 
then ask that the new function 

e i— > A afi {q e ) + G(e, b) - G(e, a) 

has zero derivative at e = 0, for all a, b. This idea leads to a version of Nocthcr's 
theorem that is at the same time widely applicable and easy to state. 



3.1 Noether's theorem in terms of BH-function and space 
change 



Definition 4. We will say that there is infinitesimal BH-invariance under space 
change q e , with G as BH-function term if 

d_ 

de' 

or, in term of the Lagrangian, if 



(A a . b (q £ ) + G(e,b)-G(e,a)) =0 Vo,6, (11) 

V / £=0 



^-(L(t, q e (t), q e (t)) + d t G(e, t)) = . (12) 

Combined with ([SJ, this invariance causes the modified function 

t^F(t) + d e G(e,t)\ s=Q 

to be a constant of motion, when q(t) solves the Lagrange equations. Of course 
we may speak of finite BH-invariance in the rare instances when the func- 
tion (|11|) is constant (Example [1] in Section [5|). 

Here is the version of Noether's theorem that uses BH-invariance under space 
change. It reduces to Theorem [T] when the BH-function is trivial. 

Theorem 2 (Noether's theorem with space change and BH-function). Suppose 
that there is infinitesimally BH-invariance under space, and that q(t) = qo(t) is 
a solution to Lagrange equation. Then the function 

t h-> d 4 L(t,q{t), q(t)) ■ d E q e (t)\ £=Q + d E G(^t)\ E=0 ■ 

is constant in t. 

Proof. As in Theorem Q] let us set F(t) = dqL(t,q(t),q(t)) ■ d e q £ (t)\ £= Q. Again 
from Lemma Q] we have that 

|- (A a , b (q e ) + G(e, b) - G(e, a)) = F(b) - F(a) + d e G(0, b) - d £ G(0, a) . 

Combining this with assumption (fTTj) . we obtain that F(b) — F(a) + d e G(Q, b) — 
d e G(Q, a) = for all a, b, which means that t H> F(t) + d £ G(0, t) is constant. □ 
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The simplest meaningful example of infinitesimal BH-invariance occurs when 
the Lagrangian function L does not depend on t. We can embed any smooth 
q(t) into its time-shift family 

q e (t) ■= q(t + e) , (13) 

whose derivatives with respect to t are the same as the derivatives with respect 
to e. Let us try a BH-function G of the form G(s, t) = eg(t). The BH-invariance 
condition (fT2|) becomes 

g'(t) = - (^(?e(*), &(*))) | e=0 = -j t L(q(t),q(t)) . (14) 
By inspection we see that (p~2]) holds with the choice 

G(e,t):=-s-L(q(t),q(t)) (15) 

We can deduce that when q{t) solves the Lagrange equations, we have conser- 
vation of energy. More on this in Example [1] in Section [S] 

3.2 Integral constants of motion 

Here is a reformulation of Theorem [2] that bypasses the concept of BH-function. 

Theorem 3 (Noether's theorem with integral constant of motion). Suppose that 
t H- q(t) is a solution to the Lagrange equation. Then the following function is 
constant: 



t i 



d 4 L(t,q(t),q(t)) -d e q e (t)\ - [ ^-L(s,q £ {s),q £ ( S )) ds . (16) 

Jt oe £ -° 



Proof. Notice that A a ^ b (q £ ) = A t0i b(q s ) — A t0ia (q E ) for all a, b, e. Take the deriva- 
tive with respect to e at e — 0, using Lemma[T]and the function F of formula |T 
on the left-hand side: 



F(b) - F(a) = —At Q ,b{qe) - —MoAle) ■ 



g- £ A t0 AQe) - J £ 

A simple separation of the terms in a and in b leads to 



Fib) ~ §- £ At A<le) = F(a) - ^ £ A tQ , a {q e ) , 

which means that the expression (|16p is constant in t. 

We may also apply Theorem [2j The condition of BH-invariance under space 
change (fT2|) is very easy to realize: with the choice G(e, t) := —A t0: t(qs) the 
function (fTTjl is trivially constant. Alternatively, if we set 

G(e,t):=-sf ^L(s,q £ (s),q £ (s)) ^ds (17) 



to 



e=0 



the function (fTT]) has zero derivative at e = 0. In either case the constant 
function 1 1-> F(t) + d £ G(e, t)\ £=0 turns out to be fl6|. □ 
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Take for example a Lagrangian of the form 



C(t,q,q):=e ht L(q,q) (18) 
and choose the time-shift space change 

q e (t) :=q(t + e). (19) 

We can compute 



(^£(t, q e (t), &(*))) | e=o = ?(*)) ■ (2°) 

An integration by part turns the BH-function (fl7|) into the following: 



d 



G(e, t) := -ee ht L(q{t) 1 q{t)) + eh f e hs L(q{s), q(s))ds . (21) 

J to 

The constant of motion is computed as 

e h %L(q(t), q(t)) ■ q(t) - e M L(q(t), q(t)) + h [ e hs L(q(s), q(s))ds , (22) 

J t 

For a mechanical interpretation of this formula see Section |6] For now it will 
suffice to notice that this conserved quantity (|22|) is nontrivial and that it de- 
pends not just on the values of t,q(t),q(t) at the current time t, but also on 
the values on the whole interval [to>t]- Hence, it is not a bona fide first inte- 
gral. Another class of examples with an integral constant of motion are the 
homogeneous potentials, as we will see in Section [7] 

3.3 The total derivative condition 

There are few and precious mechanical systems for which Noether's Theorem [3] 
does yield a true first integral, in the sense of a point function of t, q(t), q(t) that 
is constant along the motions. For the study of these systems there is a unifying 
general scheme that we are going to describe next, and that will be adhered in 
Sections [51 through [TOl devoted to examples. 

Definition 5. We will say that the total derivative condition is satisfied with 
the smooth scalar function ij)(t, q) if 



(^i(t J « s (t),« e (t)))| e=o = ^(t,9(*) ) ?(*)) Vt. (23) 

Formula (|14p . that we wrote to derive conservation of energy, is equivalent 
to a total derivative condition where ip is simply L. Compare that with (|20[) . 
whose right-hand side does not look like a total derivative. 



d 
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Theorem 4 (Noether's theorem with total derivative condition). Suppose that 
the total derivative condition is satisfied with ip, and there exists a vector-valued 
function tp(t, q, q) such that 

d £ q £ (t)\ £=0 = ^(t,q(t),q(t)) Vt , (24) 

and, as usual, that t i— > q(t) is a solution to the Lagrange equation. Then the 
following function is constant if evaluated along q(t): 

(t, q, q) i ^ OgL(t, q, q) ■ <p(t, q, q) - ij)(t, q, q) . (25) 

Proof. Simply apply Theorem [5] with G(e,t) := —e ■ tp(t,q(t),q(t j). □ 

We did not give a name to condition (IM1) because it is automatically satisfied 
whenever q e (t) is a point function of e,t,q(i),q(t),q(t ± e). More generally, a 
property of all q e (t) of this form is that when e = all their successive partial 
derivatives with respect to e and t are point functions of t, q(t), q(t), q(t). . . For 
example for the time-shift family q e (t) = q(t + e) of formulas (fT5)l and (|T5)) we 
have d e q e (t)\ E=0 = q(t). Actually, we could reach the same final conclusions 
with space variations of the more special form q e (t) = q(t) + e ■ (p(t, q(t), q(t)), 
but we have found that a delay term q(t ± e) can make formulas simpler, as we 
will see in Kepler's problem in Section [H 

While Theorem © is applicable to whatever space change q 6 (t), the total 
derivative condition is a special property of some systems, for which the space 
change q e must be chosen in a clever way. Condition (|23[) should be a con- 
sequence of the particular form of L and q £ at hand, and not of the generic 
relation (|23l) that we used in Lemma [T] 



4 Time change 

Noether's original 1918 paper [T3] did not use the concept of what we call BH- 
function, but in its place had a change of independent variables. In the present 
context this means a change of time. 

Definition 6. A time change will be a smooth real function r(e,t), defined for 
the same times t as the motion q(t), and for e close to 0, with the compatibility 
condition r(0,f) = t for all t. The trivial time change is simply r(e,t) = t for 
all t. 



4.1 (BH-)invariances under space and time change 

The idea is to take integral of the Hamiltonian action not on the fixed interval 
[a, b], but on the nonconstant interval [r(e, a), r(e, b)]. 

Definition 7. We will say that there is infinitesimal invariance under space 
and time change if 



^r(e,o),T(e,6) (<7e) 



ds 



= Wa,b, (26) 

— o 
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or, in terms of the Lagrangian, 



d_ 



(l(t(s, t), q s ( T (s, t)),q £ (T(e, t)))d t r{e, t)) 



= Vi. 



(27) 



An alternative, more common definition of invariance under space and time 
change will be discussed in Subsection 14.31 below. 

To clarify why equations (|26|) and (|27|) are equivalent, start with the action 



and perform the change of variable £ = r(e,t), keeping in mind the compatibility 
condition r(0,i) = t: 



If we take the derivative under the integral sign of this expression with respect 
to e at e ~ 0, we see that it vanishes for all a,b if and only if equation (|27p 
holds. 

Theorem 5 (Noether's original theorem with space and time change). Suppose 
that there is infinitesimal invariance under space and time change, and also that 
t i— > q(t) is a solution to the Lagrange equation. Then the following function is 
constant: 



t ^ d 4 L(t, q(t), q(t)) ■ d £ q £ (t)\ £=0 + L(t, q(t), q(t))d £ T(e, t)\ £=0 . (29) 



The proof of this Theorem will be a special case of the next Theorem [51 
While space change alone (Theorem [T]) was not powerful enough to deduce 
conservation of energy for autonomous Lagrangian systems, space and time 
change together are. Consider again a Lagrangian function L that does not 
depend on t, and take the same space change q £ (t) := q(t+e) as before, introduce 
the time change 



It is trivial to verify that ^4 T (e,a),r(e,f))( ( ?e) does not depend on e. This is a 
neat instance of finite invariance with space and time change. The induced first 
integral is the energy, again. Conservation of energy is the classical prototype 
of a conservation law that can be obtained in two different ways: (1) with BH- 
function but no time change, and (2) with time change but no BH-function. 
After Theorem [7] we will be able to exhibit a host of new examples. 

It was Bessel-Hagen [2] in 1921 who added the concept invariance up to a 
divergence (here called BH- invariance), on top of Noether's invariance under 
space and time change: 




(28) 




r(t, e) :=< — £. 



(30) 
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Definition 8. We will say that there is infinitesimal BH-invariance under space 
and time change if 

^ (A T(e , a)>T(e , 6) (<fe) + G(e, b) - G(s, a)) = (31) 

for all a, b, or, in terms of the Lagrangian, 

^ (l(t{s, t),q s (T(e, t)), q e (r(e, t)))d t r(e, t) + d t G(e, t)) = Vt . (32) 

The equivalence of equations (f3"Tj) and (|3"2"|) is left to the reader to check. 

Theorem 6 (Full Noether's theorem with space and time change and BH-func- 
tion). Suppose that there is infinitesimal BH-invariance under space and time 
change, and that also t i— > q(t) is a solution to the Lagrange equation. Then the 
following function is constant in t: 

N(t) = d 4 L(t,q(t),q(t))-d e q E (t)U=0 + 

+ L(t,q(t),q(t))d s T(£,t)\ E = + d e G(e,t)\ e = . 

Proof. We can compute the derivative in formula (|31[) using differentiation under 
the integral sign and the compatibility condition r(0,t) = t: 



= ■%-([ °' ) L(t,q £ (t),q e (t))dt + G(e,b)-G(s,a) 

0e K Jt(em) 



E=0 



= ^-! L(t,q e (t),q £ (t))dt + L(b,q(b),q(b))deT(Q,b)- (34) 

- L(a, q(a),q(a))d e r(0, a) + d £ G{0, b) - d e G(0, a) . 

To the integral term ([34|) we can apply Lemma [T] Collecting the terms in b and 
in a we obtain = N(b) — N(a) for all a, b. This means that N is constant. □ 

We leave it to the reader to find an generalization of the Definition [5] of 
the total derivative condition that is appropriate to Theorem [SJ as is done by 
Desloge and Karch [7] . In Section [TU] we will show a conservation law that is 
most simply treated with BH-invariance with both space and time change, all 
of them nontrivial. 



4.2 Equivalence of the invariances 

Infinitesimal BH-invariance under both space and time change clearly inglobates 
the other invariances, where either the gauge or the time change are absent (or 
trivial: r(e,t) = t, G = 0). After Bessel-Hagen's 1921 work [2], it apparently 
took until Boyer [5] in 1967 to notice that, given a system whose first integrals 
can be deduced with space and time change and BH-function, those same in- 
tegrals can be obtained with space change and a modified BH-function. Our 



12 



contribution in the next Theorem [7] is that those first integrals can also be de- 
rived with space change and a modified time change, with no BH-function, at 
least in our ODE setting. In Sections [5] onward the only systems to which we 
cannot apply the Equivalence Theorem [7] are the few where space change is 
enough, as in Theorem Q] 

Theorem 7 (Equivalence of invariance conditions) . Let q e (t) be a space change, 
r{e,t) a time change, G(e,t) a BH-function. Define the additional BH-function 
Q and time change T : 

Q{e,t) :=e-L(t,q(t),q(t))(der(e,t)\ 6=0 ), (35) 

T[s,t) :=g- . . (36) 

L(t,q(t),q(t)) 

(If needed, we will restrict the time t to an interval where the denominator ofT 
does not vanish). Then the following three conditions are equivalent: 



1. the infinitesimal invariance with space and time change of formula (|2T[) 
holds with t replaced by the time change r + T, i.e.: 

d I 

^{L^ qe (0,qe(0)\ S = Tie , t)+ne , t ) x 

x (d t T{e,t)+d t T(e,t))) = 0. (37) 

/ e=0 

2. the infinitesimal BH-invariance with space change of formula (| 12|) holds 
with G replaced by G + Q , i.e., 

^ (L(t, q £ (t), q e (t)) + d t G(e, t) + d t G{e, i)) |^ = 0; (38) 

3. the infinitesimal BH-invariance with both space and time change of for- 
mula (132[) holds for the given time change r and BH-function G, i.e., 

^-(L(^q £ (0,qe(0) , d t T(e,t)+d t G(e,t)) =0. (39) 



Proof. The left-hand sides of equations ()37|) . (p8| and (|39|) are identically the 
same. In terms of C(e, £) := L(d,q e (£),q c ((l)), their common value can be 
expanded out as: 

£(0, t)dl t r(0, t) + d t C{0, t)d e r(0, t) + d £ C(0, t) + dl t G(0, t) . 

This is a straightforward brute-force computation using basic two- variable chain 
rule calculus, with a little care due to nesting, and using the simplification rules 
r(0,t) = t, 9 t r(0,i) = 1. □ 



Observation 1 . Formula (|36[) is not the only possible auxiliary time change that 
makes the theorem work. An alternative choice is T(e,t) := G(e, t)/L(t, qo(t), 
qo(t)), which is possibly nonlinear in e. 
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Observation 2. One can check that the value of the constant of motion in equa- 
tion (I33p does not change if we perform either the replacements r — > r + T , 
G — > of Condition [TJ or the replacements t — >• i, G ^> H + Q of Condition [5J 

Observation 3. We may feel uneasy that formula (|36[) for T contains the La- 
grangian L at the denominator. What happens if the Lagrangian vanishes for 
some values of t? Are those values of any intrinsic importance in the trivialized 
BH-function approach? Luckily the answer is negative, thanks to this simple 
trick: choose a constant k so that L(t,qo(t),qo(t)) + k does not vanish in a 
compact interval we are interested in, and define the modified functions: 

L k = L + k, r k {e,t) := e—j- , , ■ 

L{t,q {t),q (t)) + k 

Then we may substitute the following Condition 4 for Condition[T]in Theorem[7l 



4. the infinitesimal BH-invariance with time change of formula (|27[) holds 
with L replaced by L k , the time change r replaced by r + 7fc, i.e.: 



d_ 

7k 



0. 



Of course the Lagrangians L and Lfc have the same Lagrange equations. Also, 
the value of the constant of motion in equation ([29| does not change if L is 
replaced by and r by r + 7fe- 

4.3 An alternative approach to time change 

Invariance under space and time change can be introduced in a different way. 
Let us take a time change that we write as 9 £ (t), so that we have a handy 
notation for its inverse Q~ x with respect to t. Let q £ be the composition of 
1 1-> tfe(t) with 

^(O^&o^O. g g (0 = l c °ES - (40) 

^o6» e (£) 

Compared to Definitions [7] and [51 the idea here is to take the action with variable 
cndpoints but on the reparameterized q e , instead of the original q e . 

Definition 9. We will say that there is infinitesimal alternative BH-invariance 
under space and time change if 

^ (A ee(a)Mb) {q e ) + G(e, b) - G(e, a)) | = (41) 



for all a, b, or, in terms of the Lagrangian, 



^ (i(6» e (t), q s (t), <?e(*) + 9tG(e, t) 



= Vt. (42) 

6 = 
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Formula (|4"!2)) can be deduced from (PfTj) using the change of variable £, = 9 £ (t) 
to get an integral with fixed cndpoints: 



.4 



e e (a),e e (b) 



e E (b) 



(q e ) = / L(6g e (0,$ e (0K 



S e (a) 
6 

L 



Condition (1421) is different from (|32|) . The two approaches to invariance can 
be reconciled with suitable modification of the space changes. 

Proposition 1. Suppose that a triple of space and time change and BH-function 
q £ (t), r(e, t), G{e, t) satisfies the invariance of Definition^ Define the new space 
change Qi, £ (t) in either of the following two ways: 

Q h£ {t) := q e (t - e ■ (d e r(e,t)\ e=0 )) , 
Q h£ (t) :=<fe(t)-e-(MM)le=o)g(*) 

and simply set 9 £ (t) :— r(e,t). Then the triple Qi t£ (t),9 £ (t),G(e,t) satisfies the 
alternative invariance of Definition [9l 

Conversely, suppose that a triple of space and time change and BH-function 
q £ (t),9 £ (t),G(e,t) satisfies the alternative invariance of Definition^ Define the 
new space change C>2. £ (t) in either of the following two ways: 

Q2, £ (t) :=q s (t + e-(d e r(e,t)\ e= o)) , 
Q 2 , e (t) :=<fe(t) + e-(0 e r(M)|e=o)g(t) 

and simply set r(e,i) := 9 £ (t). Then the triple Q2, £ {t),T(e,t),G(e,t) satisfies 
the invariance of Definition 



Proof. Simply notice that 

d 



de 
de 







A T ( £ ,a),T( £ ,b)(qe) _ n — 7^A? e (a),0 e (b)(Ql,e) 



6 = 







where Q 



1,6 



(a),e e (b)(<7e) 

i E o 0- 1 , and q £ := q £ o 9^ 



£ =0 ~ -Q- e A T(^a),T(sA){Q2, £ ) 



6=0 



6=0 



□ 



Theorem 8 (Full Noether's theorem with alternative space and time change 
and BH-function) . Suppose that there is infinitesimal alternative BH-invariance 
under space and time change according to Definition^ and that also t t— > q(t) 
is a solution to the Lagrange equation. Then the following function is constant 
in t: 



N(t) = d 4 L(t,q(t),q(t)) ' l*=o + (d £ 9 £ (t)\ £ ^)q(t)j + 

+ L(t, q(t), q(t))d s e e (t)\e=a + 9 s G(e, f)| e=0 . 



(43) 
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Proof. Simply apply Theorem [5] to the triple Q2, £ (t),T(e,t),G(e,t) defined in 
the second half of Proposition [TJ □ 

In terms of the alternative invariances with Definition [9l the equivalence 
Theorem [7] can be rephrased this way: 

Theorem 9 (Equivalence of invariance conditions in terms of £ ). Let q £ (t) be a 
space change, £ {t) a time change, G{e,t) a BH-function. Define the additional 
BH-function Q and time change T : 

G(e, t):=e- L(t, q(t),q(t)) (fc0 B (t)| B=o ) , T(e, t) := e ■ S^rfe: • 

L{t,q{t),q{t)) 

Then the following three conditions are equivalent: 

1. the infinitesimal invariance with space and time change of formula (|4ip 
holds with G replaced by the constant 0, the time change 9 £ (t) replaced by 
9 e {t) + T(e, t), and q £ (t) replaced by either of the following 

Qe(t) :=<?*(* + £• (9eTM)|e=0)), 

Q e {t) ~qe(t)+e-{d e T(e,t)\ e = )q(t) 

2. the infinitesimal BH-invariance with space change of formula (|12[) holds 
with G replaced by G + Q and q £ (t) replaced by either of the following 

Qe{t) :=q £ (t + E- (d e 9 e (t)\e=o)) , 

Q £ {t) :=&(t)+e - (0 e e (i)|e=o)g(t) 

3. the infinitesimal alternative BH-invariance with both space and time change 
of formula (|41l) holds for the given space change q £ , time change £ (t) and 
BH-function G. 

You will notice that, in this formulation, trivializing either time change or 
BH-function involves modifying the space change too. Also, formula (1431) for 
the constant of motion for the alternative invariance is more complicated than 
formula (|33[) . The very definition of alternative invariance involves an inverse 
time change Oj 1 , and as such is somewhat less elementary than the original 
Definition [S] These are the reasons why in this paper we have decided to stan- 
dardize time change to Definitions [7] and |S] throughout the rest of this paper, 
although Definition [5] seems to be the standard in the literature. 



5 The simplest examples 

Example 1 . Let us see how Theorem [7] works out for Conservation of Energy 
when the Lagrangian L{q,q) is autonomous. As already noted in Section [21 
there is infinitesimal BH-invariance with space change type with the choices 

q £ (t) := q{t + e) , n(e,t) := t , G x (e,t) := -eL(q{t),q{t)) ■ (44) 



1G 



Formulas ([3"5"j) and (j3"u) then become 

r n tv ^ 9eG(0,t) -L(q(t),q(t)) 

01=0, Ti(e,t):=e =g v y = -e , 

L(q{t),q(t)) L{q(t),q(t)) 

Theorem[7]says that there is infinitesimal invariance with space and time change, 
that is, with the alternative choices 

q s (t) :=q(t + e), r 2 (s,t) := n(e,t) +Ti(e,t) = t - e , G 2 = . (45) 

We have rediscovered precisely the time change t — e of formula (j30|) in Subsec- 
tion gm 

If we had started out with these last choices (j4"5"j) , we would get 

G 2 (e,t) := eL(q{t),q(t))(d E T 2 (e,t)\ E=0 ) = -eL(q(t), q(t)) , T 2 := , 

and we would be led back to infinitesimal BH-invariance with space change 
type with the original choice (fT5|) or (|44| . Whichever the approach, in the 
end the energy first integral of equation ([33)) becomes dqL(q,q) ■ q — L(q,q). 
It is unpleasant that the invariance with BH-function given by f|44[) is merely 
infinitesimal, while the invariance with nontrivial time change given by (1451) is 
actually a finite invariance. We propose here the following alternative choice 
for G 



(46) 



/t+e 

which recovers a perfect finite invariance. In fact, the quantity 

f a ,b{e) := A T2(£)a ) iT . 2 ( E>6 )(g e ) + G 2 (e,b) - G 2 (e,a) = (47) 

l-b 

L(q(( + e), g(£ + e)) d£ + G 2 (e, b) - G 2 (e, a) = 

» 

L(q(t),q(t))dt, (48) 



does not depend on e. Or again, if your favourite mnemonic reference is the 
non-integrated formula (|39[) . the expression 

L(q £ (0,qe(0) , d t n(e,t) + d t G 2 (E,t)=L(q(t),q(t)) 

does not depend on s. The associated first integral is again the energy. The 
BH-function choice (|46|) seems to be unusual, because it is an integral functional, 
and not a point function of q(t), q(t). 

We could apply again Theorem[7]to the triple q £ ,r 2 , G 2 and get a new variant 
with trivial BH-function and nontrivial time change, but we would relinquish 
the finitcness of the invariance. 

The reader can try the alternative space change q E (t) := q(t) + e • q(t) , which 
is simply the first-order expansion of q(t + e) with respect to e. 
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Example 2. (Momentum). Suppose that L(t,q,q) is invariant in the direction 
of u e K": L(t, q + eu, q) = L(t, q, q) for all J,££l, 1". Then for any 

q(t) there is obvious finite invariance for the translated family q e (t) := q(t) +eu, 
r(e, t) = t, G = 0. The constant of motion when q(t) solves the Lagrange 
equations is the component of the momentum in the direction of u: 

d 4 L(t,q(t),q(t))-u. 

Of course here Q = T = 0, and all three conditions of Theorem [7] collapse into 
one. 

Example 3. (Angular momentum) Consider a point in the R 2 plane that is driven 
by a (possibly time-dependent) central force field: L(t,q,q) := |||g|| 2 — U(t, \\q\\). 
Given a smooth trajectory q{t) in R 2 , define the rotation family 

qS) ■■= ( C0S£ ~ a[n£ )l(t), 
H v ; ^sme cose J w 

with trivial r(e,t) :— t and G = 0. It is clear that we have finite invariance: 
L(t, q £ (t), q £ (t)) does not depend on e. Noether's theorem gives the first integral 
of angular momentum for all Lagrange motions 

dqL ■ d £ q e \ £=() = q ■ ( ^ q = det(q, q). 

Again G = T=0. 

6 Dissipative systems 

Let us take up again the "dissipative" Lagrangian £(t,q,q) :— e ht L(q,q) of 
formula ([T8| . with the space change and BH- function given by (fT9| and (|21|) 
and the unusual constant of motion 

d q C(t, q(t), q(t)) ■ q(t) - C(q(t),q(t)) +h f C(s, q( S ),q(s))ds . 

For example, take m, k > 0, h = k/m, and the familiar Lagrangian L(q,q) = 
Jjto||<7|| 2 — U(q). The Lagrange equations for C are mq + S7U(q) = —kq, where 
we see a viscous resistance term — kq, which justifies our use of the term "dissi- 
pative" . The "usual" energy for the system is of course 

E{q,q) = \rn\\q\\ 2 + U(q) = d 4 L{q,q) ■ q- L(q,q) , 

and it is not conserved. Still, if we switch our attention to a "dissipative energy" 
defined as 

£(t, q, q) = e kt ' m E{q, q) = d 4 £(t, q, q) ■ q(t) - £(t, q, q) , 



18 



we can rewrite the constant of motion as 

£(t,q{t),q{t))+h f C(s,q(s),q(s))ds. (49) 

That is, the dissipative energy decreases (assuming L > 0) over time and its 
total loss over any time interval is proportional to the dissipative Hamiltonian 
action in that interval. 

A point-function first integral can be recovered in a special case as follows. 
Again with the familiar Lagrangian L(q, q) = im||<7|| 2 — U{q) we can compute 

j t (me kt/m q ■ q) = e kt/m (kq ■ q + m\\q\\ 2 - q ■ VU(q) - kq ■ q) = 



- p kt/r. 



(m\\q\\ 2 -q-VU(q)) 



Assume that U is homogeneous of degree 2, i.e., U(\q) = X 2 U (q) if A > 0. Then 
q ■ VU (q) = 2U(q), so that integrand in formula (|4"9")l becomes a total derivative: 

j t {me kt ' m q-q) = e kt ' m (m\\q\\ 2 - 2U{q)) = 2C(t, q(t), q(t)) 

and the conserved quantity (|49l) is 

e kt / m (E(q(t),q(t)) + kq(t) ■ q(t)) - e kt °/ m kq(t ) ■ q(t Q ). 

The last term containing t , as it is an obvious constant of motion, and we can 
ignore it. We finally obtain a true point- function first integral: 

e kt ' m {E{q,q) + kq-q). (50) 

Let us see how we can arrive at this same first integral with a trivial BH- 
function, following Theorem[7J The original G of formula (|21[) in our specialized 
assumptions is 

G(M) := -ee ht L(q(t),q(t)) + eh f e hs L(q(s),q(s))ds = 

J to 

k /"* 

= -eC(t,q(t),q{t))+e- / C(s, q(s), q(s))ds = 
m Jt 

= eC{t, q(t),q(t)) + e^e kt ' m q(t) ■ q(t) - » /m q(t Q ) ■ q(t ) . 
Again we ignore the to term and take the simplified 

dM) := -eC{t,q(t),q(t)) + S^ m q(t) ■ q(t) . 
The additional time change of formula (ffijf is 

-£ + le kt / m g ■ g q.q 
T{e, t):=E ^ = s + £ k -^j-j- 
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Theorem [7] tells us that we can obtain the first integral (f5U|) using 

q E (t):=q(t + e), r(e,t) := t - e + ek^ , G 2 {e,t) = 0. (51) 

Li 

The one-dimensional case, where U(q) — cq 2 , is treated by Logan [11] using 
a different set of time and space change: 

Qe(t):= (l-^jq(t-e), r 2 (e, t):=t + e , G 2 (e,t) = 0, (52) 

which he deduces by solving the generalized Killing equations. You will notice 
that the time change in ([ST]) depends on q, while the one in (|52p does not. The 
reader can compute how the equivalence Theorem [7] applies to the triple (|52[) . 



7 Homogeneous potentials 

Consider a Lagrangian of the form L(t,q,q) := im||g|| 2 — U(q), where the 
potential U is positively homogeneous of degree a, i.e., U(Xq) = X a U(q) for 
all A > 0. Define the space change as 

q e (t) = e"q(e^ a ' 2 - 1 h) . 

We can compute 

JU(t, &(*),&(*)) 



= aL + 






s=0 





Hence we can take the following BH-function 

G(s,t) := - + (f + l) ^Ld- 

and deduce the following constant of motion: 

(d 4 L ■ d £ q £ + 9 e G(e, t)) | £=0 = mq ■ q + (| - l)tE - (| + l) jf L , (53) 

where := im||(7|| 2 + U(q) is the energy, which is conserved too, because the 
system is autonomous. As you can see, the constant of motion (|53[) involves the 
Hamiltonian action and the energy. 

The special case a = —2 is specially interesting, because we get a very simple 
time-dependent first integral in the usual sense: 

F = mq ■ q — 2tE . (54) 
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Because of Theorem we can obtain this first integral with the following space 
and time change and null BH-function: 

2tL 

q s (t) = e e q{e~ 2 H), r(e,t) :=t + T(e,t) :=t + e— = (l + 2e)t. 

An alternative choice for time change is Ti(£, t) := e 2e t, which matches better 
with the space change. 

Some consequences of the first integral are easy to draw. We can take 
the antiderivative with respect to time in the equation = mq ■ q — 2tE — F and 
obtain one more time-dependent constant of motion 

F 1 = -m\\q\\ 2 -t 2 E-tF. (55) 

This means that, for any given orbit, the couple (t, \\q\\) always lie on a conic 
in R 2 : on an ellipse if E < 0, on a hyperbole if E > 0, and on a straight line 
if E = and F ^ 0, and on a single point if E = F = 0. We can also solve 
for ||# 

IMI = -Vt 2 E + tF + Fx, (56) 
m 

This formula gives exactly how the distance from the origin depends on time, 
even though we don't know the shape of the orbit. For example, when E < the 
orbit is born and dies at the origin at the instants (-F ± \JF 2 — 4EFi)/(2E). 

If we solve equation (|55|) for t and replace into equation (IMI) we obtain a 
constant of motion that does not involve the current time t, but only q(t), q(t) 
and the initial data to,q(to),q(to)- 

The central potential U(q) — —k/\\q\\ 2 is a well-known special case, for which 
Danby [6] derives formula (|55[) . Using the angular momentum, the shapes of 
the orbits can be calculated explicitly and are called Cotes' spirals [HI p. 83]. 

8 Kepler's problem 

In all the examples that we have seen so far, the total derivative condition holds 
without requiring the Lagrange equation on q(t). In this Section we will see 
examples where we need the Lagrange equation, or at least some consequences 
of them. 

The Lagrangian function of Kepler's problem with its associated Lagrange 
equation are 

where k > is a parameter. Notice that for all Kepler motions q is parallel to q. 
We will assume that q is a vector in M. 2 to simplify some formulas. 

As for all autonomous Lagrangians, the energy is conserved. The Lagrangian 
is also invariant under rotations around the origin. Noether's theorem gives the 
first integral of angular momentum det(g, q). 
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Kepler's system enjoys also a recondite infinitesimal invariance. Let us take 
again a smooth trajectory q(t) in R 2 \ {(0,0)} (not necessarily a solution to 
Lagrange equations, for now), a vector u £ M 2 , and define the family 

q e (t) := q(t) + (q(t) ■ u)q(t + e) - (q(t + e) ■ u)q{t) = 
= q(t)+det(q(t),q(t + e))u 1 - 

with u 1 - — (J "q 1 )u. It is clear that qo(t) = q(t). This family q e is different and 
simpler than the one found in the literature (for example, see Levy-Leblond [10] , 
formula (36)), whose formula would be 

q(t) + |(2( g (i) • u)q{t) - (q{t) ■ u)q(t) - (q(t) ■ q(t))u 

The reader can check that the following relations hold: 

^L(q e (t), q e (t)) | s=Q - det(q(t), q(t)) det(tt, q(t)) + -| ( fc ^yj) • 
If we define the BH-function 

G(e,t):=-ekA^, (59) 

and if q(t) is any motion for which q(t) is parallel to q(t), then we have infinites- 
imal invariance of the form (1121): 



d_ 

ds 



L(q e (t),q e (t))+d t G(e,t) 



= 0, 

e=0 



which is covered by Theorem [21 condition [3j with the trivial r(e,t) = t. Also, 
q(t) and q(t) are parallel whenever q(t) is a solution to Kepler's equation ([57]) . 
Therefore Noether's Theorem [5] (or Th. 0} yields the following constant of mo- 
tion: 

{q-u)\\q\\ 2 - {q-u)(q- q)- k——. 

h\\ 

Since the vector u £ M 2 is arbitrary, we have the vector- valued first integral 

QUI] 2 - (Q ■ q)q ~ k rr-rr , 

\\q\\ 

which is called the Laplace-Runge-Lenz vector. 

According to Theorem [71 if we define the additional time change T as 

r(M);=e m . u 



L(q(t),q(t)) \\q(t)\\L(q(t),q(t)) ' 

there is infinitesimal invariance also under the nontrivial time change (e, t) i— > 
t + T{e, t), without BH-function. 
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Let us show a simple example where Noether's theorem is applied to a partic- 
ular solution, leading to a nontrivial function that is constant along that single 
solution, but not along most others. Consider the following family of uniform 
circular motions in the plane with the same period but a phase shift: 

q e (t) := e ue R(cos(u)(e + t)), $m(oj(e + t))) , 

where R > 0, u — ^k/R 3 . The function t i-> q e (t) is a Kepler motion only 
when e = 0, as we check at once. Still, there is infinitesimal invariance with 
trivial r(e, t) = 0, G = 0: 

9 r i , . . , x \ d i R 2 uj 2 ._ k \ 

If we apply Noether's theorem we obtain that the (square of) the speed t i-> 
1 1 qo (£)|| 2 is constant along the circular motion. The speed is clearly a nontrivial 
function that is not constant along any Kepler motions except circular ones. 

9 New superintegrable systems 

A recent paper 19J introduced the following Lagrangian systems in two dimen- 
sions: 

L(x,y,i,y) = xy - g(x)y, 

and exhibited some classes of the function g for which there is either weak 
Lyapunov instability or an isochronous center. What matters here is that those 
isochronous cases exhibit the rare property of being superintegrable, because 
they have three independent hrst integrals, out of four degrees of freedom. Here 
we are going to use Noether's theorem framework to compute new g classes that 
lead to super-integrability. 

Two first integrals are obvious, and they do not need any assumption on g: 
^x 2 + V(x) and yx + g(x)y, where V is any primitive of g. Notice that the 
Lagrange equations are 

x = ~g(x), y = -g'(x)y. (60) 

Let us search for a third constant of motion starting from the following 
space-change family: 

q e (t) = (x(t) + ef(x{t),x(t)),y{t)y 

where / is a function to be determined. We can compute, using also the La- 
grange equations (|6U|) : 

^L{qS)Ae{t))\ s=Q =yf-g'(x)yf = j t {yf) - (f + g'(x)f)y. 
If we impose that 

f + g'(x)f = (61) 
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then we achieve the total derivative condition, there is infinitesimal BH-invari- 
ance with the BH-function G(e,t) — —ey(i)f, and Noether's theorem gives the 
following first integral 

d 4 L(q(t),q(t))-deq s (t)\ s =o + d e G(e,t)\ e = = yf-yf. (62) 

It is the first case so far where the BH-function G is a function of q, and not 
just of t and q. One more example will be in the next Section. 

The condition (|61[) is a partial differential equation in the unknown func- 
tion fix, x). We will search for smooth solutions that are even functions of x: 

f(x, x) = a{x) + x 2 P{x) + x^hix, x 2 ) (63) 

for new unknown functions a{x), /3(x), h(x, x). Let us replace the ansatz (|63p 
into the equation (fBTj) and set x = 0. We get 

= a(x)g'(x)-g(x)a'(x) + 2g(x) 2 p(x) = g(x) 2 {2^x) + -^-(a(x)/g(x))) (64) 

This suggests to introduce the new function fi(x) = a{x) / g{x), so that equa- 
tion becomes simply /? = —fJ.'/2. We can change the ansatz into the 
more special form 

f(x, x) — g(x)fi(x) — x 2 /j,'(x) + x A h(x, x 2 ) . 
Let us try with the choice h = a constant. Equation (|6"Tj) becomes 

iiV» - 6ag'(x)) + l -x 2 (2Aag{x) 2 + 2»(x)g"(x) + 

+ 3g'(x)(i'(x)-3g(x)ti"(x)) =0, 

which splits into the equations 
H"'(x) - 6ag'(x) = 

2Aag{x) 2 + 2fi(x)g"{x) + 3g'{x)n'(x) - 3g(x)n"(x) = , 

to which we can apply the standard local existence results for ODEs. Any 
solution of this equations will give a function g for which the system is superin- 
tegrable. The additional first integral (1521 becomes 

a(4:X 3 yg(x) + x A y) - ^x 3 yn"(x) + ^x 2 yfi'(x) - xy(i{x)g'(x) + yg(x)n{x) . 

It can be checked that the special case a = 0, with the initial data g(Q) = 0, 

g'(0) = 1, gives exactly the class of systems that were found in Section 5, 

entitled "Explicit superintegrable systems", of the original paper [19 , for which 

the additional first integral is a polynomial of degree 3 with respect to x. The 

case o^O brings superintegrable systems with a first integral containing a term 
• 4 

m x . 

With some patience it is possible to pursue these ideas further, for example 
by studying the case when h(x, x 2 ) is of the form j(x) + ax 2 . 
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10 Particle in a plane-wavelike external field 



The following time-dependent Lagrangian is taken from a paper [4 by Bobillo- 
Ares: 

L(t,q,q) = ^\\q\\ 2 -U{q-tu), q,q£R n , 

where u is a fixed vector in K™ and U a smooth potential. The associated 
Lagrange equation is 

q + VU(q-tu)=0 (65) 

In terms of the energy 

E(t, q, q) = 8gL{t, q,q)-q- L(t, q, q) = ^\\q\\ 2 + U(q - tu) , 

it is easy to check that q ■ u — E is a first integral. Let us see how we can deduce 
it from Noether's theorem in our framework. Starting from a smooth q(t) and 
following Bobillo-Ares, we introduce the following space and time changes: 



q E (t) = q(i) +su, r(e,t):=t + e 
Let us try infinitesimal invariance: 
d 



(66) 



d t r(e,t) 

§- £ (\W + e)f - U(q(t + e)+eu-(t + e)u)) 
q(t + e) ■ q\t + s) - VU(q(t + e) - tu) ■ q(t + e) = 



(t + e)- ^2q(t + e) - (q(t + e) + VU(q(t + s) - tufj^j = 
\\q(t + e)\\ 2 - q(t + e) ■ (q(t + e) + VU(q(t + e) - tu)} . 



_ d_ 

~ dt 

When e — this expression becomes 
d 



j t \\q(t)\\ 2 -q(t)-(q(t) + \7U(q(t)-tu) 



which further reduces to 

d_ 

de 



\ f=r(e,t) / e=0 Ctt 



if q(t) solves Lagrange equations (|65|). This means that we have infinitesimal 
invariance as in Theorem [6] with the following choice of BH- function: 

G(e,t) := -e\\q(t)\\ 2 , 

and the first integral ([33]) given by Noether's theorem is 

q(t)-U + L(t,q(t),q(t))-\\q{t)\\ 2 , 
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as expected. We do not know another example where it is natural enough to 
have infinitesimal invariance with space change, time change and BH-function, 
all of them nontrivial. Also, notice that this is one more example where the BH- 
function G depends on q, so that is not of the more familiar form s • tp(t, q(t)). 

According to Theorem [71 if we define the additional time change and BH- 
function 

Q{e, t):=e- L(t, q(t),q(t)) (9 £ r(e, i)| e=0 ) = e ■ L(t, q(t), q(t)) , 

Tie t) = z 9eGM = e Uit)r 
{,) - L(t,q(t),q(t)) L(*,g(t), ?(<))' 

we can attain infinitesimal invariance also with either of the alternative choices 

rx(s,t) := r[e,t) +T(s,t) , Gi{e,t) = 0, 
T 2 (s,t) :=t, G 2 (e,t) :=G{e,t)+Q(e,t), 

for the same space change q e . 

A possible alternative choice for the BH-function term in (|66|) is the following 

-\\q(t + e) - q(t)\\ 2 /e if e ^ 
if e = , 

which is not linear in e, and is not a point function of q(t),q(t). 



G 3 (s,t) 
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